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SECTION  I:  GENERAL  PLASMA  THEORY  AND  SIMULATION 


A. 


Vortex  Formation  and  Particle  Transport 
in  a  Cross-Field  Plasma  Sheath 


K.  Theilhaber  and  C.K.  Birdsall 
Electronics  Research  Laboratory, ' 
University  of  California, 
Berkeley,  CA  94720 


Abstract 


The  time-dependent  behavior  of  a  transversely  magnetized,  two-dimensional 
plasma-wall  sheath  has  been  studied  through  particle  simulations,  which 
have  shown  that  the  cross-field  sheath  develops  into  a  turbulent  boundary 
layer,  with  strong  potential  fluctuations  and  anomalous  particle  transport. 
The  driving  mechanism  is  the  Kelvin-Helmholtz  instability,  which  arises 
from  the  sheared  particle  drifts  created  near  the  wall.  The  sheath  acquires 
an  equilibrium  thickness  Ij  ~  5  p,-,  and  maintains  large,  long-lived  vortices, 
with  amplitudes  64>  ~  — 2r,/c,  which  drift  parallel  to  the  wall  at  roughly 
half  the  ion  thermal  velocity.  The  sheath  also  maintains  a  large,  spatially- 
averaged  potential  drop  from  the  wall  to  the  plasma,  in  sharp  distinction 
with  the  unmagnetized  sheath,  where  the  plasma  potential  is  higher  than  at 
the  wall.  Accompanying  the  vortices  is  a  spectrum  of  shorter-wavelength 
fluctuations  which  induce  an  anomalous  cross-field  transport.  A  central 
simulation  result  is  that  for  Wp,  >  2u>ci,  the  transport  scales  like  Bohm 
diffusion,  a  result  for  which  we  have  a  qualitative  analytic  model. 


Submitted  to  Physical  Review  Letters,  June  1988. 
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B.  Theoretical  Studies  of  Edge  Physics  and  PlasmaAVall  Interactions* 

C.K.  Birdsall  and  W.S.  Uwson  UCB 

D.Hamed  SAIC  San  Diego 

The  proposed  woik  represents  a  collaboration  between  SAIC  and  the  University  of  California,  Berkeley 
to  address  important  issues  in  the  physics  of  the  plasma  edge.  In  particular,  state  of  the  art  tools  will 
be  applied  to  understanding  the  coupling  between  the  plasma  edge,  modeled  by  a  fully  kinetic  particle 
code  and  the  plasma  center  modeled  by  a  semi  implicit  resistive  MHD  code.  This  work  is  expected  to 
significantly  enhance  our  understanding  of  divertors  and  the  plasma  H  mode.  In  addition,  the  manage¬ 
ment  plan  for  this  collaboration  includes  significant  interactions  with  experimental  groups  involved  in 
these  issues.  (Excerpt  from  letter  firom  Dr.  D.A.  Hitchcock,  Fusion  Thetxy  BraiKh,  EXDE  September  14, 
1987.) 

The  original  plan  was  that  collaboration  between  fluid  and  particle  researchers,  SAIC  and  UCB 
respectively,  could  lead  to  substantial  advances  in  understanding  magnetized  plasmas  near  walls,  such 
as  found  in  tokamaks,  with  limiters  and  diverters.  After  the  original  plan,  interest  at  DOE  increased  in 
learning  more  about  edge-plasma  behavior  at  and  near  an  ICRF  antenna  with  electrostatic  shield,  as 
used  for  tokamak  heating;  we  added  this  objective,  upgrading  the  complexity. 

The  Berkeley  effort,  with  a  small  doctoral  student  group,  was  originally  slated  to  begin  May  1, 
1987  which  became  about  October  1,  1987.  (Support  was  not  received  in  Berkeley  until  months  later.) 
We  had  the  good  fortune  that  Dr.  William  Lawson  agreed  to  stay  around  a  bit  after  his  thesis  was  com¬ 
pleted,  in  order  to  start  some  of  the  plasma-edge  work  proposed.  The  work  performed  so  far,  directly 
on  this  contract,  is  mostly  by  him.  Unfortunately,  having  Just  begun,  we  were  then  told  we  would  only 
receive  six  months  support,  to  end  April  1,  1988. 

ACCOMPLISHMENTS 

Our  basic  machine  for  doing  our  computer  experiments  on  a  magnetized  two-dimensional  plasma 
in  contact  with  a  wall  is  our  many-particle  code  ES2  developed  here  by  Dr.  Kim  Theilhaber.  (We  are 
also  developing  PDW2,  using  much  of  ES2,  along  the  lines  of  our  well  known  and  widely  used  Id 
code,  PDWl.)  The  major  recent  work  with  this  code,  by  Dr.  Theilhaber,  is  "Vortex  Formation  and  Par¬ 
ticle  Transport  in  a  Cross-Field  Plasma  Sheath"  Memorandum  No.  UCB/ERL  M88/21,  20  March  1988, 
100  pages,  widely  circulated.  A  short  version  .has  been  submitted  to  Physical  Review  Letters  and  a 
more  detailed  version  is  to  be  submitted  to  the  Physics  of  Fluids. 

*  This  is  a  progress  report  only.  A  final  report  is  being  prepared  with  all  of  the  detailed  results. 
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This  report  goes  into  great  detail  on  the  sheath  formation,  edge  waves,  instability  (Kelvin- 
Helmholtz),  vortex  formation  and  coalescence,  and  anomalous  Bohm-Iike  transport  of  particles  to  the 
wall.  This  effort,  which  took  place  over  two  years,  with  joint  DOE  and  other  support,  covers  the 
undriven  magnetized  plasma  sheath  very  well.  It  forms  one  of  the  starting  points  for  the  joint  SAIC- 
UCB  Berkeley  work,  which  will  now  be  described. 

It  is  to  be  noted  most  firmly  that  while  fluid  codes  and  implicit  particle  codes  may  be  applied  to  a 
magnetized  plasma  near  a  wall,  such  codes  are  not  cabable  of  resolving  all  of  the  detailed  physics 
needed,  on  the  scale  length  of  a  Debye  length  or  a  gyroradius,  in  the  sheath.  The  latter  physics,  which 
is  needed  for  the  transport  across  the  sheath  (particles,  energy,  heat)  may  only  be  produced  with  any 
amount  of  trust  by  explicit  particle  codes  such  as  we  are  using.  Such  physics  is  needed  for  these  fluxes, 
as  well  as  the  wave  and  instability  descriptions.  It  is  the  marriage  of  particle  codes  near  the  wall  and 
fluid  codes  in  the  bulk  plasma,  including  the  use  if  implicit  coding  away  fi’om  the  wail,  that  is  pro¬ 
posed. 

UNDRIVEN  PLASMAAVALL  RESULTS 

In  the  plasma-wall  region  which  we  will  label  the  sheath  ,  we  find  an  unusual  plasma,  with: 

n,  *  n, 

gradient  lengths  scaling  with  or  p, 
f(v)  non-Maxwellian 
momentum  reversals  (of  one  species) 
surface  modes  (e.g.,  Kelvin-Helmholtz  waves) 
wall  charges  and  currents 

wall  absorption,  reflection,  emission  (i.e.,  neutrals) 


e(|> 

KT, 
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This  un-plasma  is  very  different  from  the  usual  nearly  uniform,  nearly  quiet,  nearly  Maxwellian,  nearly 
neutral  textbook  plasma;  the  un-plasma  tends  to  demand  kinetic  treatment,  using  particles. 

Four  long  computer  runs  were  made  on  the  ES2  codes  of  Theilhaber,  starting  with  an  appreciable 
vacuum  gap  between  the  metal  absorbing  wall  and  the  plasma  edge.  The  object  was  to  observe  whether 
there  was  any  appreciable  new  physics  with  the  plasma  initially  away  from  the  wall  or  up  against  the 
wall.  The  latter  case  was  studied  in  detail  by  K.  Theilhaber,  leading  to  Kelvin-Helmholtz  instability, 
formation  of  vortices,  coalescence,  and  Bohm-like  diffusion  to  the  wall. 

The  conclusions  are  thac 

(a)  the  plasma  in  all  runs  moves  across  the  vacuum  gap  of  many  p, ,  perhaps  Sp, ; 

(b)  the  inidal  instability  is  a  finite  gain  ion  gyioradius  drift  wave,  developing  into 

what  may  be  a  Kelvin-Helmholf^  velocity-shear-driven  mode. 

Hence,  as  was  an  objective  of  these  runs,  it  appears  that  the  plasma  tends  to  fill  the  space  alloted. 
Therefore,  it  is  quite  all  right  to  start  studies  with  the  plasma  against  the  wall,  to  be  described  now. 

THE  ANTENNA-DRIVEN  PLASMA-WALL  RESULTS 

In  the  previous  section,  we  found  that  the  plasma  would  rather  readily  close  up  a  gap  of  a  few 
ion  gyroradii,  suggesting  that  it  is  permissible  to  start  the  plasma  right  up  against  the  wall.  Hence,  we 
felt  free  to  start  the  driven  plasma  runs  without  the  vacuum  gap.  Long  runs  were  made,  with  similar 
parameters  as  before. 

We  found  that  the  electrostatic  near  field  of  the  ICRF  antenna  substantially  increases  the  particle 
flux  and  energy  to  the  antenna  sheild,  which,  in  the  laboratory,  results  in  a  large  increase  in  sheild 
sputtering.  It  is  this  increase  (due  to  the  above  and  other  effects  as  well)  which  limits  operation  in 
some  of  the  large  tokamaks;  such  effects  are  well  stated  in  the  summary,  by  RJE.  Aamodt,  of  the 
ICRF/Edge  Physics  Workshop,  Boulder  CO,  March  30  —  April  1,  1988. 

Obviously,  many  more  runs  and  some  more  diagnostics  are  called  for,  as  well  as  invention  in  the 
electrostatic  shield  and  antenna  configuration,  so  as  to  alter  (reduce)  the  wall  sputtering. 


_4_ 


SECTION  n:  PLASMA  WALL  PHYSICS,  THEORY  AND  SIMULATION 


A.  A  Proposed  Particle-In- Cell  Method  for  Modeling  Small 
Angle  Coulomb  Collisions  in  Plasmas 


S.  E.  Parker 

Electronics  Research.  Laboratory 
Univ.  of  California,  Berkeley 

December  6,  1988 


1  Introduction 

We  are  proposing  a  computational  method  to  self-consistently  model  small  angle  coUisional 
effects.  This  method  may  be  used  in  standard  Particle-In-Cell  (PIC)  plasma  simulations 
to  include  collisions,  or  as  an  alternative  to  solving  the  Fokker-Planck  (FP)  equation  using 
finite  difference  methods.  The  distribution  function  is  represented  by  a  large  number  of 
particles.  The  particle’s  velocities  change  as  a  function  of  the  drag  force,  and  the  diffusion 
in  velocity  is  represented  by  a  Wiener  process.  This  is  similar  to  previous  Monte-Carlo 
methods  [1,2],  e.xcept  we  calculate  the  drag  force  and  diffusion  tensor  self-consistently.  The 
particles  are  weighted  to  a  grid  in  velocity  space  and  the  associated  “Poisson’s  equations” 
are  solved  for  the  Rosenbluth  potentials.  First  the  approximation  for  small  angle  Coulomb 
collisions  is  discussed.  iNe.xt,  the  FP-PIC  collision  method  is  outlined.  Then  we  show  a  test 
of  the  particle  advance  modeling  an  electron  beam  scattering  off  a  fixed  ion  background. 
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2  Small  angle  Coulomb  collisions 


The  FP  equation  for  describing  small  angle  Coulomb  collisions  can  be  solved  numerically 
using  finite  difference  techniques.  An  alternate  method  is  to  follow  the  evolution  of  a  large 
number  of  particles.  Using  the  notation  of  Trubnikov  [3],  an  infinitestimal  “cloud”  of  test 
particles  can  be  represented  by  the  quantities  <  An,-  >  and  <  Au.Avj  >  to  the  same  order 
of  accuracy  as  the  FP  equation[3].  These  are  defined  as: 


<  Au,  >= 


dvi(t) 

dt 


(1) 


<  AviAvj  >= 


(2) 


The  bars  represent  averages.  These  may  be  interpreted  as  ensemble  averages  over  many 
initial  states,  or  as  a  straight  average  by  letting  the  number  of  particles  in  a  cloud,  N,  get 
large,  N  oo.  We  follow  the  second  interpretation  allowing  a  simpler  simulation  method. 
Then  r,  is  the  average  velocity  of  the  local  cloud.  The  i  and  j  subscripts  are  the  3  velocity 
components  in  cartesian  coordinates,  (ui  =  Vx.,V2  =  Uy,  t>3  =  r*).  <  Au,  >  is  the  drag  force 
felt  on  an  average  particle  in  the  cloud.  <  Au,  At>j  >  is  the  spreading  or  diffusion  in  velocity 
space  of  the  cloud. 

Numerically  we  change  the  velocity  of  the  particles  by  using  the  following  equation  [4]: 


Auj(t")  =<  Au,  >  Af  +  5,jAWi(t") 


(3) 


where  B,j  satisfies: 


BikBjk  =<  AviAvj  > 


(4) 


AW,(t'’)  =  )  and  W,(t)  is  a  vector  function  that  represents  a  Wiener  process 

(or  Brownian  motion)  having  the  following  properties  <  W,(f)  >=  0  and  <  W,(t)Wj(f)  >  = 
tSij.  B  and  W  are  not  unique  and  only  need  to  satisfy  the  above  criteria.  In  our  model  we 
choose: 

AW, it’')  =  sfAtR,  (5) 


where  R,  are  independent  random  numbers  satisfying  <  it,  >=  0,  and  <  R,Rj  >=  6,j.  We 
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have  found  a  a  matrix  B  such  that  BB"^  =  A,  given  by: 


( 


Bij  = 


0 

0 


(6) 


fe-tefe 

By  letting  A,j  =<  Au.Awj  >  we  note  that  equation  (6)  satisfies  eq.  (4)  and  can  therefore 
use  this  B  in  equation  (3). 

The  “field”  quantities  are  obtained  from  two  functions  (j>  and  ip,  the  Rosenbluth  poten¬ 
tials,  which  in  turn  are  solved  by  the  two  Poisson’s  Equations  [3]: 

VV  =  f(vi)  (7) 

V  V  =  <i>^{vi)  (8) 

where  /3  is  the  superscript  representing  the  field  species,  /?  =  (t,  e)*  Then  <  An,  >  and 
<  AviAvj  >  are  obtained  in  terms  of  the  Rosenbluth  potentials  using  the  following  equation 
(see  reference  [l]  for  a  derivation): 


<  An, 


i/0  ^^0 

dVi 


(9) 

(10) 


=  X  ( — j  is  a  constant  (given  here  in  cgs  units),  and  A  is  the  Coulomb  logarithm. 
\  / 

Q  and  0  represent  the  test  and  field  particles,  respectively  (e.g.  a  =  (i,e)  and  0  =  (i,e)). 


3  An  Outline  of  the  Method 

For  the  particle  advance  we  start  with  the  simplest  scheme  [1,2]  similar  to  Euler’s  method 
except  there  is  an  added  diffusion  term  (the  last  term  on  the  right): 

<  Av.  >"  At  -I-  B^jVAtRj  (11) 

Where  Rj  are  independent  random  numbers  having  the  two  properties  given  above.  Follow¬ 
ing  the  methodology  of  PIC  simulation,  we  weight  the  particles  to  a  grid,  calculate  the  field 
quantities,  advance  the  particles  and  repeat  the  process.  The  basic  algorithm  is  as  follows: 
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Step  1.  Weight  the  paxticles  to  a  grid  in  velocity  space  using  linear  interpolation. 
Step  2.  Solve  equations  (6)  and  (7)  on  the  grid  for  <i>  and 

Step  3.  Obtain  <  Auj  >"  and  B^-  on  the  grid  using  equations  (6),  (9)  and  (10). 

Step  4.  For  each  particle,  obtain  <  Av,  >"  and  by  interpolating  from  the  grid  to 
the  particle  location  u,-.  Then  update  the  velocity  using  equation  (11). 


4  Solution  to  the  “Field”  Equations 


In  step  2.  above,  we  need  to  solve  Poisson’s  equation  twice  in  3  dimensions.  That  is, 
=  3,  where  u  =  and  s  =  (/,<A).  To  simplify  the  problem,  we  assume  azimuthal 

symmetry  in  and  work  in  a  2  dimensional  cylindrical  geometry  (r,z)  where  r*  = 
and  z  =  u,.  Since  the  system  is  independent  of  the  spatial  coordinates  we  define:  x  =  v*, 
y  =  Vy  and  z  =  v,  without  confusion.  A  finite  Fourrier  transform  (FT)  is  made  in  z,  then 
we  can  solve  for  each  k  using  the  following  finite  difference  scheme  [5]: 


ArSrj 


where: 


rj+i/2  =  rj  +  -Ar 
=  ’’1+1/2  ~  ’’>-1/2 

K  is  the  finite  difference  analog  to  fc[5]: 

"  =  (I;)  (^) 


(12) 


(13) 

(14) 


(15) 


This  tridiagonal  system  is  solved  for  each  k,  and  an  inverse  FT  is  made  to  get  u,  u{r,z)  — 
FT-’(ii(r,jb)). 

The  next  step  is  to  calculate  ..  Avi  >  and  <  AviAvj  >  from  (p  and  tp.  Using  the  change 
of  coordinates  and  the  assumption  of  azimuthal  symmetry  we  obtain  the  following; 


d4> 

dv. 


{—<l>r,-<l>r-,4>z 

r  T 


) 


(16) 
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f  +  (i  -  ^)A  ^Ar  -  ^A  fA. 

^Ar-^A  ^Ar  +  i7-^)A  ^Az  (17) 

jAx  ^Az  Az  ) 

We  approximate  A^A^A,  At,  and  Az,  by  central  differencing  <j>  and  V*.  Then  use  equations 
(8)  and  (9)  to  get  <  Au,  >  and  <  Au.Auj  >  on  the  grid. 


5  Test  of  the  particle  advance 

As  an  initial  test  of  the  particle  advance,  eq.  (11),  we  model  an  electron  beam  scattering 
off  infinitely  massive  cold  ions.  For  this  test  case  <  Av,-  >  and  <  Av,  Auj  >  are  given  by 
the  two  simple  analytic  expressions: 

<  At;.-  >=  -C^  (18) 


where  C  =  - .  We  do  not  calculate  the  field  quantities  on  the  grid,  but  rather,  use 

Air 

equations  (18)  and  (19)  explicitly.  This  gives  a  good  test  of  the  particle  advamce,  equation 
(11).  We  expect  the  total  x- momentum  to  be  given  by: 


Px, totality  —  Px,total{^  —  0)€Xp  ^ 


where  vq  is  the  initial  beam  velocity.  For  the  run  shown:  C  =  1,  vq  =  1,  vxe  =  0.01  and 
At  =  0.001.  Figure  1  shows  the  initial  beam.  Figure  2.  is  the  total  x-momentum  vs.  time. 
Figures  3  and  5  show  snapshots  of  (v^,  Vy).  Figure  4  and  6  also  are  (v^,  Vy)  scatter  plots  but, 
a  small  slice  in  Uj,  jujl  <  0.1.  The  energy  error  is  less  than  0.2  percent  and  the  momentum 
error  is  less  then  1  percent. 

6  Future  Work 


We  are  developing  a  code  to  test  the  FP-PIC  method  outlined  above.  In  the  future  we 
plan  to  implement  the  algorithm  in  a  Id  3v  electrostatic  code  to  study  combined  collective 
electrostatic  and  collisional  effects.  In  this  case  we  would  assume  spatial  homogeneity  for 
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calculating  the  coUisional  terms.  If  this  assumption  is  not  valid  one  would  have  to  partition 
f{v)  into  spatial  regions  j,  and  calculate  fj{v)  for  each  spatial  region.  One  would  have 
to  ensure  that  enough  particles  in  each  spatial  zone  to  adequately  fill  out  the  distribution 
function. 

In  addition,  we  need  to  address  the  accuracy  and  stability  issues  associated  with  this 
method,  and  study  the  feasibility  of  using  a  more  accurate  particle  advancing  scheme  than 
that  given  by  equation  (10). 
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Figure  1:  (urt^v)  beam  at  t=0. 
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Figure  2;  x-momentum  vs.  time. 


-11- 


vx 

Figure  3:  (u^Wy)  for  the  beam  at  t=0.1. 
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vx 

Figure  4:  Slice  in  u,  of  beam  at  t=0.1,  |t;j|  <  0.1. 
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Figure  6:  Slice  in  of  beam  at  t=l,  |u.|  <  0.1. 
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B.  Direct-Implicit  Particle  Simulations  of 

Collisionless  and  Collisional  Bounded  Plasmas. 

R.  J.  Procassini  and  C.  K.  Birds  all 
and 

B.  I.  Cohen  (Lawrence  Livermore  National  Laboratory) 

Direct-implicit  particle-in-ceU  simulation  methods  are  applied  to  the  study  of 
botmded  plasmas  systems,  with  the  goal  of  obtsuning  a  self-consistent  description 
of  plasma-surface  interactions  in  the  region  near  the  divertor  plate  of  a  tolmmak. 
Several  important  phenomena  in  this  region  of  the  scrape-ofF-layer  (ionization  of 
recycled  neutral  gas,  chaxged/neutral  collisions,  etc.  )  have  characteristic  time  and 
space  scales  that  are  long  compared  to  the  electron  plasma  frequency  and  Debye 
length.  Direct-implicit  methods  provide  an  efficient  means  for  accurately  reproduc¬ 
ing  these  effects,  without  having  to  resolve  electron  relevamt  time  and  space  scales. 
We  will  present  the  results  of  large  Up^At,  Ax/X^g  implicit  particle  simulations  of 
a  bounded  plasma  system  in  both  the  coUisionless  emd  collisional  limits.  These 
simulations  are  not  capable  of  resolving  A^e,  hence  the  sheath  width  is  based  upon 
the  gild  spacing.  The  total  potential  drop  through  the  presheath  and  sheath  in  this 
implicit  limit  will  be  investigated,  and  compeired  to  theoretical  predictions. 
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c. 


Self-Consistent  Calculations  of  the  Sheath 
and  Presheath  in  a  Bounded,  Collisionless 
Plasma  with  Different  Source  Distributions. 


R.  J.  Procassini  and  C.  K.  Birdsall 
and 

B.  I.  Cohen  ( Lawrence  Livermore  National  Laboratory) 

Two  different  models  for  the  calculation  of  the  sheath  and  presheath  in  a  bounded, 
coUisionless  plasma  have  recently  been  proposed  by  Emmert,  et.  al.  '  and  Bissell 
and  Johnson^.  These  models  differ  in  the  choice  of  the  source  distribution  func¬ 
tion  for  the  ions.  Emmert  uses  fs{v)  =  \v\f\f{v),  while  Bissell  and  Johnson  use 
fs(v)  =  where  fM(v)  is  a  Maxwellian.  Both  models  assume  Boltzmann  elec¬ 

trons.  We  employ  particle-in-cell  techniques  to  study  the  effect  of  these  different 
source  distribution  functions  on  the  sheath  and  presheath,  without  making  the  as¬ 
sumption  of  Boltzmauin  electrons.  Our  electron  distribution  function  is  therefore 
cut-off  by  the  collector  sheath  potential.  Self-consistent  potentieJ  profiles  for  each 
type  of  source  distributions  will  be  presented. 

[1]  G.  A.  Emmert,  et.  al.  ,  Phys.  Fluids^  23,  803  (1980). 

[2]  R.  C.  Bissell  and  P.  C.  Johnson,  Phys.  Fluids,  30,  779  (1987). 
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A. 


SECTION  m;  CODE  DEVELOPMENT 


Particle- Circuit  Simulation  of  Traveling-Wave- Tubes 

Ian  J.  Morey 

A  code  is  being  developed  to  simulate  a  traveling- wave- tube.  It  uses  the  particle-in-cell 
technique  to  model  the  electron  beam,  combined  with  a  leapfrog  scheme  for  modeling  the 
coupled  transmission  line.  An  advantage  of  the  particle-in-cell  techmque  is  the  ease  with 
which  space-charge  effects  can  be  modelled.  Since  the  code  was  initially  intended  to  be 
used  as  a  teaching  aid,  it  is  being  developed  to  run  on  a  PC  with  many  diagnostics  of 
the  interaction  displayed  graphically.  However,  it  is  believed  that  the  code  could  become 
a  useful  research  tool  if  run  on  a  faster  computer.  At  present,  msmy  parameters  can  be 
changed  while  the  code  is  running,  such  as:  source,  line  and  load  impedances,  coupling 
between  the  beam  and  the  transmission  line,  space-charge  effects  and  the  input  signal 
amplitude.  Other  parameters  can  only  be  changes  at  the  stzirt  of  each  simulation,  such  as: 
beam  energy,  beam  current,  tube  length  and  wave  velocity. 


B. 


ESI  on  IBM-Compatible  Personal  Computers 
John  P.  Verboncoeur 

The  Electrostatic  1  Dimensional  Periodic  Plasma  Simulation  (ESI)  code  originated  by  A. 
Bruce  Langdon'  in  FORTRAN  on  the  MFE  Cray  computers  and  subsequently  translated  into  C 
by  T.  Lasinski  has  been  ported  to  the  MS-DOS  environment.  The  PC  version  of  ESI  adds 
several  features  to  the  Cray  version  which  are  described  below. 

The  MS-DOS  version  of  ESI  can  simulate  periodic  plasmas  comprised  of  up  to  8  different 
species.  The  total  number  of  particles  for  unmagnetized  simulation  is  16000,  and  up  to  8000 
for  magnetized  species.  PC  ESI  has  the  capability  to  weight  charge  to  up  to  8192  grid  points 
using  zero  order  (nearest  grid  point)  or  1st  order  (particle-in-cell)  weighting  schemes.  The 
code  can  run  up  to  1024  timesteps,  storing  electrostatic,  kinetic,  total,  and  Fourier-decomposed 
mode  energies.  All  the  standard  input  parameters  of  ESI  are  supported  (except  plot  intervals; 
PC  ES 1  plots  diagnostics  interactively  with  a  keypress).  The  input  parameters  are  contained  in 
an  input  deck  similar  to  the  namelist  input  deck  of  the  Cray  version.  PC  ESI  is  capable  of 
moving  1000  particles  per  second  on  an  IBM  AT,  2000  on  an  IBM  PS/2  Model  60,  and  40(X) 
on  a  20  MHz  80386  PC. 

Accessibility.  The  MS-DOS  version  of  ESI  runs  on  any  IBM-compatible  PC,  AT,  or  PS/2. 
The  proliferation  of  these  machines  enables  almost  any  student  or  researcher  to  run  and  modify 
ESI.  We  have  already  distributed  over  50  copies  of  the  code,  and  it  has  been  used  in 
mini-courses  and  offsite  demonstrations.  The  complete  source  code  for  ESI  is  distributed  with 
the  program  to  facilitate  additional  diagnostics,  new  particle  pushers  and  field  solvers,  etc. 

Phase-Space  Animation.  PC  ESI  displays  a  real-time  phase  space  animation  with  each  species 
in  a  different  color.  Using  the  animation,  one  can  observe  one  dimensional  collisions, 
accelerations,  and  particle  trapping.  PC  ESI  can  also  leave  traces  along  the  trajectories  of  the 
particles  in  phase  space. 


1  C.  K.  Birdsall  and  A.  B.  Langdon,  Plasma  Physics  Via  Computer  Simulation  (1985) 


-17- 


c. 


PDWl  on  IBM  Compatible  Personal  Computers 
John  P.  Verboncoeur 

The  Plasma  Device  Workshop  One  Dimensional  Bounded  Plasma  Simulation  (PDWl)  is  being 
ported  to  the  MS-DOS  environment.  The  original  code,  due  to  William  S.  Lawson^,  was 
written  in  FORTRAN  for  the  MFE  Cray  computers.  In  1986,  Tom  L.  Crystal  translated  the 
FORTRAN  source  code  into  C.  We  now  intend  to  create  a  version  of  PDWl  on  the  PC 
analogous  to  PC  ESI.  The  initial  version  will  be  completed  in  the  4th  quarter  of  calendar  1988. 


2  William  S.  Lawson,  PDWl  User’s  Manual,  Memorandum  UCB/ERL  M84/37  (1984). 
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D. 


ESI  for  Windows 

John  P.  Verboncoeur  and  Vahid  Vahedi 

Sinailar  to  PC  ESI,  Windows  ESI  adds  a  number  of  features.  Using  the  multiple  display 
windows  capabilities  of  Windows,  Windows  ESI  can  display  any  or  all  of  the  diagnostics 
concurrently  with  the  phase  space  animation.  The  diagnostics  are  updated  with  each  step  of  the 
animation,  resulting  in  real-time  display  of  space-charge  density,  fields,  etc.  The  windows  can 
be  sized  and  moved  to  Suit  the  user’s  needs.  The  mouse -driven  graphical  interface  is  well 
suited  to  the  beginner  or  casual  user. 

In  addition,  Windows  ESI  adds  pop  up  menus  for  diagnostics  as  well  as  a  menu-driven  help 
system.  We  expect  to  submit  Windows  ESI  in  the  4th  quarter  of  calendar  1988,  satisfying  the 
requirements  of  the  IBM  Windows  Software  Development  Contract 


-19- 


E.  The  Construction  of  Distribution  Functions 


Communication  from  Dr.  Andrew  Zachary*  to  BJ.  Cohen  LLNL 


Your  notes  on  the  construction  of  the  “true”  distribution  function  are  much 
appreciated.  As  we  have  already  discussed,  one  possible  method  of  construction 
/  requires  uniform  binning  in  emonuniform  variable.  For  example,  to  determine 
f{v±)v±dvx,  the  bins  are  distributed  uniformly  in  ux-  Similarly,  for  f(v)v^dv, 
the  bins  are  uniformly  distributed  in  v^.  Your  notes  and  the  different  figures  in 
TESS  show  that  you  have  adopted  this  approach.  However,  nommiform  variables 
reqiiire  an  unacceptable  numb^  of  bins  whenever  the  particles  have  a  large  range 
in  velocity.  In  my  cosmic  ray  simulations,  the  particles  have  momenta  which  lie  in 
the  intend  (0, 200),  so  I  cannot  use  yoxir  technique. 

There  is  an  alternative  -  bin  the  particles  uniformly  in  a  tmiform  variable,  and 
then  try  to  “xindo”  the  binning.  For  convenience,  let  me  use  v  to  denote  the  velocity 
variable,  and  suppose  I  have  a  series  of  uniformly  spaced  velocity  bins  t  =  1,  . . 
N.  For  each  point  Ui,  I  have  an  integral  equation  for  /  of  the  form 


=  / 


Sid(v  —  vi)f(v)v‘‘dv. 


Here,  F{vi)  denotes  the  value  returned  by  the  usual  distribution  function  routine 
as  the  estimate  of  the  distribution  function  at  v,-.  Also,  5id(u  -  Vi)  is  the  usual 
one-dimensional  linear  weight  function 


Sid{v  -Vi)  = 


Vi-i<V<Vi 

-  W)  Vi  <v  <  u,+  i 


which  represents  the  shape  of  an  individual  particle.  Furthermore,  v°  is  the  velocity 
space  volume  factor,  where 

{0  if  u  =  v„  U|| 

1  if  u  =  Ux 

2  if  V  =  Ivl 


•  The  University  of  Chicago 
Department  of  Astronomy  and  Astrophysics 
5640  South  Ellis  Avenue 
Chicago,  Illinois,  60637 
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As  it  stands,  Eq.  1  is  intractable  without  some  method  of  approximating  /(v). 
After  some  experimentation,  I  believe  that  the  “best”  solution  involves  approximat¬ 
ing  /  on  the  interval  by  a  piecewise  linear  function 


/(v.-i)  +  -^{fM  -  v.-i  <v<Vi 

/(w.+i)  -  ^(/(*'.)  -  /(^‘+i))  Vi  <v<  U,+i 


This  approximation  places  no  constraints  on  the  function  form  of  /,  although  it 
does  require  that  df/dv  >  <Pf/dv^  in  each  subinterval  (u,-, u.+i).  Except  for  possi¬ 
bly  pathological  distribution  functions,  the  piecewise  linear  approximation  seems  a 
reasonable  dioice. 

Now  the  integral  in  Eq.  1  can  be  done  in  closed  form.  After  some  tedious 
algebra,  the  result  is  a  tridiagonal  set  of  coupled  equations 

F{vi)  =  -f-  Ci^if{vi)  Ci^i+if{vi+i) 


where 


The  symbol  6a, 2  is  the  usuzd  Kronecker  delta. 

At  the  upper  and  lower  boundaries,  the  coeflBcients  C,,,  are  modified  as 


Au 

6 


gVlAV  -b 


12 


Av  f  2,1 
3  2 


ujAv  -b 


'  ^ 

6 

1  A  . 

Cn,n  =  gUyvAu  -I - 


(M 

10 


a  =  0 
a  =  1 
0  =  2 

0  =  0 
0  =  1 
0  =  2 


My  test  calculations  show  that  the  resulting  set  of  tridiagonal  equations  are 
stable  and  robust  over  a  wide  range  of  velocities.  Furthermore,  the  reconstructed 
distribution  function  seems  very  close  to  the  expected  “true”  distribution. 
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I  have  sxifficient  confidence  in  this  approach  that  I  have  extended  it  to  handle 
the  reconstruction  of  a  2D  distribution  function.  Siwpoae  I  need  to  reconstruct 
/(up,w±)  or  /(A»»y)  fro™  momenta  F(t;||,t;x)  or  F(/i,u).  Again,  assume  that  F 
is  known  on  a  set  of  uniformly  spaced  mesh  points  {xi,yj),  with  i  =  I,  . . N  and 
j  =  I,  M.  The  2D  eqmvalent  of  Eq.  1  is 

Vj)  =  /  f  S2d(x  -  *i.y  -  Vi) fix,  y)x^y*dzdy.  (2) 

Jxi-l 

I  am  going  to  simplify  the  problem  by  awniTiiing  that  x  is  “unweighted”,  i.e.  that 
0  =  0.  Fortunately,  —  Xi,y  —  yj)  factors  into  two  parts 

S^ioix  -Xi,y-  yj)  =  Sid(x  -  Xi)SiDiy  -  yj). 

Then  the  inteCTal  in  Eq.  2  couples  the  value  of  /  at  (xj,  yj)  to  the  values  of  /  at  the 
8  nearest  neighboring  points. 


-  *,;  +  !  — 

*  +  j  + 1 

© 

*  1  1 

*  -^1  -/ 

t  -f  1,; 

© 

i-lj-l  - 

1 

1  ' 

1 

In  the  spirit  of  the  one-dimensional  technique,  I  approximate  /(z,  y)  by  a  piece- 
wise  bilinear  function.  Explicitly,  in  each  region  shown  above,  I  write 

Wfi{x,y)  =  (Xi-t-i  -x)(y_,+i  -  y)/(ii,yj)  +  (xj+i  -  x)(y  -  yj)/(z„y,+i) 

+  (x  -  Xi)(y_,+i  -  y)/(x,+i,y,)  -H  (x  -  z,)(y  -  y,)/(z,+i,y^+i) 

W^/2(r,y)  =  (x  -x._i)(y,+i  -y)/(x,-,y,)-|-(x-Xi_i)(y-y,)/(zi,yj+i) 

+  (li  -  x){yj+i  -  y)/(x,_i,y,)  +  (xi  -  x)(y  -  y>)/(x,_i, y_,+i) 

^f3ix,y)  =  (x  -  Xi_i)(y  -  y>-i)/(xi,y,)-f  (x  -  x,_i)(y_,  -  y)/(x,-, yy_i) 

+  ixi  -  x)(y  -  y_,_i)/(z<_i,y,)  -f-(xi  -  x){yj  -  y)/(ii_i, yj._i) 

^/4(x,y)  =  (xi+i  -x)(y-y,_i)/(zi,y,)-l-(xi+i  -  x)(yj  -  y)/(xi,  yj_i) 

+  ix-  Xi)(y  -  yy_i)/(ii+i,y>)-f-(x  -  Zi)(yy  -  y)/(z,+i,yy_i) 

where  W  =  Ax  Ay.  After  some  more  tedious  algebra,  the  integrals  in  each  region 
can  again  be  done  in  closed  form.  The  results  are 

Q  =  1 
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(yy  +  +  2^^^  ■*■  ^XAj+i  +  2-^i+i.i+i) 

^2(x,-,  yj)  =  —  (yy  +  -^)(/i,j  +  2^i-ij)  + 

■F3(x<,y>)  =  y  (yy  -  ^)(/i,y  +  “  ^KAi-i  +  ^/i-i.i-i) 

i^4(x„ yy )  =  y  [(yy  -  ^)(/.-.y  +  i/i+,,y)  +  ^(yy  -  ^)(/i,y-i  +  \fi-x,j-l) 

and 

a  =  2 

^’.(x.v.)  =  ^[(v?  +  ^  (/<.i  + 

+  2  Vj^y  +  "  jQ~^  (/>.j+l  +  2/*+i.>+i)] 

f2(ii,vi)  =  y  [(y?  +  ^  (/i,i  +  |/i-ij) 

+  2  (y?  +  y^^y  +  )  (/ij+i  + 

+  5  (y/  -  y/Ay  +  (/fj-i  +  \u-i.i-i) 

F.(x.,  „,)  =  E  [(,j  _  yyA,  ^  ^ 

+  2  ^yj ~  vi^y  +  ^ iQ^  'j  ^ 

When  the  four  terms  are  added  together,  the  resiilt  is 

■^(xi,yy)  =  ■^ij/i-lj-hl  +  ^i.y/i.y+l  +  ^«>7/»+l,ji+l 

+  •^•.y/i+i.y  "h  ^i,jfi—ij—i  +  ^i,y/i,y— 1  ■!■  ^«,y/i+i,y— i 

Now,  if  I  let  £  =  (j  —  +  t,  then  I  can  tremsfonn  these  equations  a  single, 

giant  matrix  equation.  This  matrix  has  9  non-zero  diagonals,  and  eacK  row  has 
the  pattern:  3  non-zero  elements,  followed  by  N  —  1  zeros,  followed  by  3  non-zero 
elements,  then  another  N  —  1  zeros,  3  more  non-zero  elements,  and  then  the  rest  all 
zeros.  I  have  not  yet  written  the  routines  to  solve  the  resulting  sparse  set  of  linp-ytr 
equations.  I  am  talking  with  some  of  the  math  people  at  MFE,  and  hope  to  have 
the  routines  ready  by  the  end  of  the  month. 
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SECTION  IV:  SHEATH  SEMINAR 


^  University  of  Californio,  Berkeley  CA  94720  5  May  1968 

FINAL  SCHEDULE  for  seminar  on 

PLASMA  SHEATHS  AND 
OTHER  BOUNDARV  LAYERS 


This  is  on  interdisciplinary  plasma  seminar,  sponsored  by  faculty  from  the 
departments  of  Chemical  Engineering,  Electrical  Engineering  and  Computer 
Science,  Mechanical  Engineering,  Nuclear  Engineering,  and  Physics.  The 
organizer  is  Prof.  Charles  K.  Birdsall,  EECS,  191M  Cory  Hall,  (415)643-6631. 
See  the  initial  flyer  for  more  description.  Partial  support  for  this  seminar 
comes  from  a  MICRO  project  with  Varian  gift,  gratefully  acknowledged. 

Meetings  are  Mondays,  2-4pm,  January  25  through  May  9,  in  240  Bechtel  Engr 
(No  meetings  on  February  15  or  March  26  due 'to  University  holidays.) 

The  order  being  followed  is  roughly  1/3  DC  sheaths,  1/3  RF  driven  plasmas, 
and  1/3  magnetized  plasmas  and  2d  effects.  A  major  objective  is  to  make 
ties  among  theories,  experiments,  and  simulations  for  various  models,  the 
references  will  have  labels  T,E,S. 

January  25  Prof  C.K.  Birdsall,  EECS,  UCB,  review  of  the  clossicl  theory  and 
experiments  on  the  plasma  potential  of  L.  Tonks,  I.  Langmuir, 

T,E  "A  General  Theory  of  the  Plasma  of  an  Arc”,  Physical  Review  34, 
September  15,  1929,  pp.  876-922. 

February  1  Dr.  Lou  Ann  Schwager,  Sandia  (Livermore)  and  UCB,  a  review  of 

the  works  of  the  1960‘s  and  1970’s  leading  up  to  her  own  theory 
and  simulations  (dissertation,  1987), 

T,E  "Kinetic  Model  of  the  Collisionless  Plasma  Region,  Applied  to 

Thermionic  Convertors,  Q-Mechines,  and  the  Boundary  Layer  of  a 
Maxwellian  Plasma  (covering  both  source  and  collector  sheaths)" 

February  6  Julia  Little,  Math  UCB  and  Scott  Parker,  NE  UCB  will  review 
further  theory  and  computation,  building  on  the  two  previous 
talks,  where  more  complete  solutions  have  been  computed,  from 
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the  following  references: 

T  S.  A.  Self,  "Exact  Solution  of  the  Collisionless  Plasma  Sheath 
Equation",  Physics  of  Fluids  6,  December  1963,  pp.  1762-1768. 

T  S.  A.  Self,  "Asymptotic  Plasma  and  Sheath  Representations  for 

Low-Pressure  Discharges"  Jour.  Appl.  Phys.  36,  Feb.  1965,456-459 

T,E  D.  Dunn  and  S.  A.  Self,  Static  Theory  of  Density  and  Potential 

Distribution  in  a  Beam-Generated  Plasma"  Jour.  Appl.  Phys.  35, 
January  1964,  pp.  113-122. 

T  J.  Parker,  "Collisionless  Plasma  Sheath  in  Cylindrical  Geometry" 
Phys.  Fluids.  6,  1963,  pp  1657-1658. 

T  G.  Emmert,  R.  Wieland,  A.  Mense,  and  J.  Davidson,  "Electric  Sheath 
and  Presheath  in  a  Collisionless,  Finite  Ion  Temperature  Plasma", 
Phys.  Fluids  23,  April  1980,  pp.  603-812. 

T,E  S.  Meassick,  M.H.  Cho,  and  N.  Hershkowitz,  "Measurement  of 

Plasma  Presheath"  IEEE  Trans,  on  Plasma  Science,  PS- 13,  April 
1985,  pp.115-119. 

February  22  Leo  Eskin,  Stanford  University  (with  Prof.  S.  A.  Self)  "Electrode 

T  Plasma  Boundary  Layer  and  Current  Transfer  for  a  Weakly- 

Ionized,  Collision  Dominated  Thermal  Plasma,  Including 
Generation" 

February  29  Blake  Wood,  EECS  Dept.,  UCB.  Review  of  boundary  conditions  on 

T  ion  energy  and  velocity  entering  collector  sheath,  starting  with 

Bohm  (1949),  continuing  through  Harrison  end  Thompson  (1959), 
Hall,  Hall  and  Bernstein,  Chodura,  Franklin,  Lamb  (1987). 

March  7  (a)  Prof.  M.  A.  Lieberman,  EECS  Dept.,  UCB.  "Analytic  Solution  for 
T  RF  Driven  Sheaths  “ 

T  (b)  Chris  Goedde,  Physics,  UCB. "  Calculations  of  RF  Sheath 
Electron  Heating" 

March  14  Dr.  Ian  Brown,  LBL.  Plasma  Immersion  lOn  Implantation  Sheaths 

March  21  Ajit  Paranjpe,  Stanford  University  (with  Prof.  S.  A.  Self),  "RF 
T  Plane  Parallel  Low-Pressure  Discharge  in  SF5,  Solving  the 

Boltzmann  Equation  for  fg(t)  and  Calculating  the  Charged-Particle 
Balance  Equations  for  Electrons,  *  and  -  Ions". 
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April  4  (a)Prof.  David  Graves,  Chem.  Engr.  UCB,  "Review  of  Modeling  for  Short 
T,  S  Mean  Free  Path  RF  Discharges",  including  works  of  Groves,  Godyok. 

(b)  Prof.  Sidney  Self,  ME  Dept.,  Stanford,  "Boundary  Conditions  for 
T  Continuum  Equations  of  the  Diffusion  Type,  with  and  without 
Migration  dut;  to  Body  Forces*  with  Applications  to  Diffusion  in 
Neutral  Gases  and  Piosmos,  Brownian  Diffusion  (Diffusiophoresis), 
Turbulent  Diffusion,  Neutron  Diffusion,  Photon  Diffusion  (Radiative 
Diffusion). 

April  1 1  (a)  George  Mislum,  EEC5  Dept,  UCB,  "The  Complete  Design 
T  Procedure  for  Symmetric  RF  Discharges" 

T,E,5  (b)  Dr.  Ian  Morey,  EECS  Dept.,  UCB,  "Short  Review  of  the  Workshop 
on  ICRF-Edge  Plasmas"  held  at  Boulder  CO,  March  30-April  l . 

April  18  (a)  Prof.  Reiner  Stenzel,  Physics  Dept.  UCLA,  "Instabilities  and 
E  Emission  Lines  of  Sheath  Plasma  Resonances" 

T  (b)  Prof.  S.  Self,  Stanford,  continuation  of  April  4  talk  on 
boundary  conditions. 

April  25  (a)  Arthur  Sato,  Physics  Dept.  UCB,  "The  Plasma  Resonance  Probe, 
T,E  with  works  of  Harp,  Kino,  Pavcovich  et  al." 

T,S  (b)Richard  Procassini,  NE  Dept.,  UCB,  The  One-Dimensional 
Magnetized  Sheath",  based  on  works  by  R.  Chodura  et  al. 

May  2  (a)  Dr.  Steve  Savas,  Applied  Materials  Corp,  Santa  Clara." 

TEM  Surface  Waves  in  13.56MHz  Plasma  Process  Discharges" 

(b)  Dr.  Kim  Theilhaber,  UCB,  now  LLNL,  "Theory  and  Simulations 
of  the  Two-Dimensional  Magnetized  Sheath,  including 
Kelvin-Helmholtz  Instability  and  Growth  to  large  "Steady" 
Vortices" 

May  5  (Special  seminar,  out  of  sequence,  visitor) 

Dr.  Rod  Boswell,  Australian  National  Univ.,  Canberra,  "Plasma 
Processing  with  Inductively  Coupled  RF  Plasma  Discharges" 

May  9  (a)  Dr.  Ian  Morey,  EECE  Dept.,  UCB,  "Particle  Simulations  of 

Parallel  Plate  RF  Discharges"  (works  done  at  ANU,  Canberra) 

(b)  Prof.  C.K.  Birdsall,  EECS  Dept,  UCB  "Review  of  Sheath  and 
Boundary  Layer  Models,  as  presented  in  this  seminar. 

(c)  End  of  the  seminar  party,  somewhere  nearby. 
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25  January  1988 


B.  PLASMA  SHEATHS  AND  OTHER  BOUNDARIES  Seminar 

CHARACTERIZATIONS 


We  will  be  reviewing  a  few  dozen  works  this  semester,  with  a  main 
objective  of  finding  which  sheath  and  boundary  layer  models  work  and  which 
do  not,  as  learned  from  theory,  experiment  and  simulation.  We  will  be  wise 
to  follow  some  kind  of  characterization  in  these  reviews,  looking  for  the 
details  of  the  models.  Hence,  some  such  will  now  be  suggested,  and 
hopefully  followed,  making  comparisons  easier. 

Beneral  characterizations: 

number  of  dimensions:  1 , 2,  or  3 
time  independent  (dc) 
time  dependent  (ac) 
not  magnetized 

magnetized  (usually  applied  dc  Bq) 
driven  (battery,  RF,  other) 

Particle  characterizations  (electrons.  Ions,  neutrals) 
sources,  sinks,  at  walls,  in  the  volume 
after  creation,  how  treated. 

electrons  particles,  fluid,  Boltzmann,  cut-off  Maxwellian,  other 
ions:  similarly,  or  cold,  warm,  beam,  other 
neutrals;  temperature,  pressure,  excited  states,  other 
equations  of  motion  include  terms  for  (list) 
equations  of  continuity  include  terms  for  (list) 

Wall  characterizations  (metal,  dielectric,  vacuum,  another  plasma,  etc.) 
floating,  no  net  currrent 

connected  to  external  circuits  and  dc/ac  sources  by  wires,  antennas, 
waveguides,  optical  fibers 
absorbing,  reflecting,  emitting 
sputtering,  etching,  depositing 

Characterization  of  results  (from  theory,  experiments,  simulations) 

potentials,  fields,  densities,  pressures,  currents,  with  critical  values 
for  stable/unstable  regions,  changes  of  state,  other 
wave  propagation,  linear,  nonlinear 

oscillations,  fluctuations,  special  features  (like  virtual  cathodes) 

fluxes  of  particles,  energy,  heal  to  walls 

scalings  obtained  (e  g.,  particle  flux  to  wall  goes  as  I/Bq),  or  in 

lengths,  lime,  amplitudes. 
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These  characterizations  are  meant  as  guidelines,  sort  of  a  recipe  for 
making  notes  while  reading  the  articles  and  preparing  transparencies  for 
presentations.  They  are  intended  to  make  the  task  efficient  as  well  as 
suggesting  a  common  mode  for  the  presentations.  Lastly,  they  give  us  a 
ready  way  of  comparing  the  models  made  and  tested  by  various  authors  and 
aid  in  pinpointing  room  for  Improvement,  Please  give  them  a  try. 

Our  colleague  5,  Kuhn  in  Innsbruck  (plus  visits  here)  published  the  chert 
below  in  one  of  his  works,  showing  a  particular  set  of  characterizations. 
You  might  find  making  such  charts  useful  also. 


TABLE  II.  Comparison  of  some  representative  contributions  to  the  kinetic 
theory  of  linear  longitudinal  plasma  oscillations  in  one'dimensional  geome¬ 
try.  _ 


Ok  O  O  —  0»  »IO 

r»  t—  A  roloo 

»  o  »  a  » 


«e  I « I  CM  I I  oo  I  n 


c. 


University  of  Califomio,  Berkeley 


14  January  1988 

This  is  a  list  of  references  begun  by  C.  K.  Birdsan*for  use  in  our  Spring 
1988  interdisciplinary  plasma  seminar, 

PLASMA  SHEATHS  AND  OTHER  BOUNDARY  LAYER  DYNAMICS 

The  references  are  given  by  year,  author,  and  journal,  with  some 
connotations,  like  theory,  experiment  and  or  simuation  Cr,E.S),  and  other 
comments.  We  will  draw  on  the  list  already  begun  by  our  Innsbruck 
colleague,  S.  Kuhn,  and  by  our  Plasma  Theory  and  Simulation  Group  members, 
especially  Lou  Ann  Schwager. 

We  are  interested  in  unusual  plasma  regions,  with  large  potential 
energies,  e<ji  >  KT,  with  appreciable  charge  imbalance,  n^  f  ng,  with 
relatively  steep  gradients  in  densities  Oq  g  ),  In  potentials  and  fields  <j>,  E, 
and  In  temperature  T.  We  expect  wall  charges,  wall  absorption,  emission  and 
reflection,  along  with  coupling  to  external  circuits,  sources.  The  plasma 
regions  may  have  elastic  and  ineleastlc  collisions  (scattering,  exciting. 
Ionizing,  radiating,  etc.)  as  well  as  generation  of  particles  by  all  means. 

Such  regions  show  up  In  sheaths,  at  probes,  in  double  layers,  at  the  edges 
of  freely  expanding  plasmas,  both  magnetized  and  not.  All  laboratory  and 
space  plasmas  have  some  kind  of  boundaries  someplace,  most  of  which  will 
Influence  the  Interior  behavior  to  some  extent.  These  regions  ore  labelled 
unusual  only  In  the  sense  that  they  do  not  obey  the  contraints  which  are 
commonly  put  on  plasma  models  In  order  to  make  the  models  tractable. 

Indeed,  it  may  be  that  what  we  are  labelling  unusual  may  be  more  common. 

It  may  well  be  that  the  final  understanding  of  plasmas  may  come  from 
understanding  what  Is  really  going  on  at  the  plasma  boundaries.  It  is 
certainly  true  that  what  we  know  about  plasmas  is  obtained  by  making 
contact  through  plasma  boundaries. 

*Note;  The  list  attached  is  not  the  "final  list"  that  I  plan  to  prepare  during 
the  seminar.  It  is  a  copy  of  works  collected  by  Dr.  Siegbert  Kuhn,  our 
associate  in  Innsbruck.  As  you  will  note,  there  are  only  5  from  the 
1920’6  through  the  ISSO's;  there  are  50  in  the  1960’s,  34  in  the  1970’s 
and,  so  far,  86  in  the  1980 ’s.  Please  bring  additions  to  my  attention. 


Kuhn’s  list  plus  Birdsall’s  sheath  collection 
are  now  typ^  in  REFER  on  our  local  computer. 
We  will  publish  this  list  soon. 
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SECTION  V:  JOURNAL  ARTICLES,  REPORTS,  TALKS,  VISITS 


Journal  Articles 


Niels  F.  Otani  and  Bruce  I.  Cohen,  "Effect  of  Large-Amplitude  Perpendicularly  Propagating  Radio  Fre¬ 
quency  Waves  on  the  Interchange  Instability,"  Phys.  Flui^,  31,  pp.  158-176,  January  1988. 

Niels  F.  Otani,  “A{^lication  of  nonlinear  dynamical  invariants  in  a  single  electromagnetic  wave  to  the 
study  of  the  Alfven-ion-cyclotron  instability,”  Phys.  Fluids,  31  June  1988,  pp.  1456-1464. 


Reports 


K.  Theilhaber,  "Vortex  Formation  and  Particle  Transport  in  a  Cross-Field  Plasma  Sheath,"  University  of 
California,  Beriteley,  Memorandum  No.  UCB/ERL  M88/21.  March  20,  1988. 

Lou  Ann  Schwager  and  C.  K.  Birdsall,  "Collector  and  Source  Sheaths  of  a  Hnite  Ion  Temperature 
Plasma,"  University  of  California,  Berkeley,  Memorandum  No.  UCB/ERL  M8/23,  April  13,  1988. 

Lou  Ann  Schwager,  "Effects  of  Secondary  Electron  Emission  on  the  Collector  and  Source  Sheaths  of  a 
Finite  Ion  Temperature  Plasma,"  University  of  California,  Berkeley,  Memorandum  No.  UCB/ERL 
M88/24,  AprU  13,  1988. 

Lou  Ann  Schwager,  "Effects  of  Ion  Reflection  on  the  Collector  and  Source  Sheaths  of  a  Finite  Ion 
Temperature  Plasma,"  University  of  California,  Berkeley,  Memorandum  No.  UCB/ERL  M88/25,  April 
13,  1988. 

Poster  Papers 

ICRF/Edge  Physics  Workshop,  March  30-April  1,  1988,  Boulder,  Colorado. 

W.  S.  Lawson  and  C.  K.  Birdsall,  "Undriven  Plasma  Wall  Interaction  Simulations,  Showing  Tur¬ 
bulence  with  and  without  an  Initial  Vacuum  Gap." 

W.  S.  Lawson  and  C.  K  Birdsall,  "Antenna  Driven  Plasma  Wall  Interaction  Simulation,  Showing 
Local  Turbulence  and  About  3  Times  Larger  Flux  to  the  Wall." 

IEEE  Conference  on  Plasma  Sciences,  June  6-8,  1988,  Seattle,  Washington.  (Abstracts  follow) 

W.  S.  Lawson,  M.  A.  Lieberman,  and  C.  K.  Birdsall,  "Electron  Dynamics  of  RF  Driven  Parallel 
Plane  Reactor." 

A.  Friedman,  S.  L.  Ray,  C.  K.  Birdsall,  and  S.  E.  Parker,  "Multi-Scale  Particle-in-Cell  Plasma 
Simulation:  Timestep  Control  Criteria  and  some  Tests." 


Talks 


C.K.  Birdsall,  "Plasma-Sheath-Surface  Dynamics  via  Particle  Simulations,"  June  9-10,  1988,  at  NSF 
Workshop  on  New  Directions  in  Plasma  Engineering. 
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is  the  driving  frequency  u,  but  there  is  a  persistent  small  oscillation 
at  Up  —  see  ^(<)  plot  below  for  i  =  600  to  800; 

The  peak  voltage  is  at  the  predicted  amplitude  when  0,  =  0,  but 
increases  slowly  with  the  thermal  velocity  so  that  at  V|  =  2  it  is  30% 
higher  than  predicted. 
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